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ABSTRACT 

Cosmological (ACDM) TreeSPH simulations of the formation and evolution of galaxy 
' groups and clusters have been performed. The simulations invoke star formation, chem- 

QQ \ ical evolution with non-instantaneous recycling, metallicity dependent radiative cool- 

, ing, strong star-burst, and (optionally) AGN, driven galactic super-winds, effects of a 

CO ' meta-galactic UV field and thermal conduction. Results for two clusters, one Virgo- 

\ like (T~3 keV) and one (sub) Coma- like (r~6 keV), are presented. At z—Q the stellar 

■ contents of both clusters consist of a central dominant (cD) galaxy surrounded by 
' cluster galaxies and intracluster stars. The intracluster (IC) stars are found to con- 
tribute 20-40% of the total cluster B-band luminosity at z=Q and to form at a mean 

Q^' redshift Zf ~3 being on average about 0.5 Gyr older than the stars in cluster galaxies. 

UBVRIJHK surface brightness profiles of the IC star populations are presented; the 
profile of the larger cluster matches the observed V-band profile of the cD in Abell 1413 
(T ~8 keV) quite well. The typical colour of the IC stellar population is B-R=1.4-1.5, 
^ \ comparable to the colour of sub-L* E and SO galaxies. The mean Iron abundance of 

the IC stars is about solar in the central part of the cluster (r ~100 kpc) decreasing 
to about half solar at the virial radius. The IC stars are a-element enhanced with a 
weak trend of [0/Fe] increasing with r and an overall [O/Fe]~0.4 dex, indicative of 

■ dominant enrichment from type II supernovae. 

\ The IC stars are kinematically significantly colder than the cluster galaxies: The 

velocity dispersions of the IC stars are in the inner parts of the clusters (r ^100-500 
kpc) only about half of those of the cluster galaxies increasing slightly outward to 
about 70% at r=l-2 Mpc. The typical projected velocity dispersion in the Virgo-like 
cluster at ii> 50 kpc is 300-600 km/s depending on orientation and projected distance 
from the cluster center. Rotation is found to be dynamically insignificant for the IC 
stars. The velocity distributions of IC stars and clusters galaxies are in one cluster 
highly radially anisotropic, in the other close to isotropic. 

Key words: cosmology: theory — cosmology: numerical simulations — galaxies: 
clusters — galaxies: formation — galaxies: evolution 
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1 INTRODUCTION 

Current hierarchical or "bottom-up" structure formation 
theories predict that the field star populations of haloes 
of galaxies like the Milky Way should consist partly of 
stars originally born in small proto-galaxies and later tidally 
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stripped from these by the main galaxy or through inter- 
action with other proto-galaxies. The fraction of halo field 
stars of such origin depends (apart from obvious cosmic vari- 
ance effects) on the properties of the dark matter (e.g., cold 
dark matter, CDM, vs. warm dark matter, WDM: for WDM 
the hierarchy is broken at the free-streaming mass scale, 
e.g., Sommer-Larsen & Dolgov 2001) and the nature of the 
star-formation process. For example, some halo stars may be 
born in rapidly expanding, supernova driven proto-galactic 
supershells and hence be detached from their parent proto- 
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galaxies right from the outset in space, as well as phase-space 
(e.g., Mori ct al. 1997; Sommer-Larsen, Gotz & Portinari 
2003). If dark matter is cold and star- formation in expanding 
supershells is not a common phenomenon the above fraction 
could bo quite large, in principle at least, approaching unity 
(e.g., Helmi et al. 2003). The halo stars resulting from tidal 
stripping or disruption of a proto-galaxy will stay localized 
in phase-space for a long period and several such "streams" 
of halo stars have been detected in the haloes of the Milky 
Way and M31 (e.g., Holmi ct al. 1999; Ferguson ct al. 2002). 

From the point of view of structure formation, galaxy 
clusters can be seen as scaled up versions of galaxies in 
hierarchical scenarios. In particular, gravity is expected to 
strip or disrupt cluster galaxies in a similar way as satellite 
galaxies in galaxy haloes and a population of intracluster (or 
"field") stars should thus reside between the cluster galax- 
ies. It has been known for decades that cD galaxies often 
are embedded in extended envelopes presumed to consist of 
stars tidally stripped off galaxies in the process of being en- 
gulfed by the cD (e.g., Oemler 1976; Dressier 1979). In recent 
years it has been possible to perform quantitative studies of 
these stellar envelopes through ultra-deep surface photom- 
etry of the general stellar population (e.g., Gonzalez ot al. 
2000; Gonzalez, Zabludoff & Zaritsky 2004; Feldmeier et al. 
2002, 2004a), or imaging/spectroscopy of individual plane- 
tary nebulae (e.g., Arnaboldi ct al. 2002, 2003; Feldmeier et 
al. 2004b) or supernovae Typo la (Gal- Yam et al. 2003). The 
potential importance of intracluster stars in relation to the 
chemical enrichment of the intracluster medium (ICM) was 
recently discussed by Zaritsky, Gonzalez & Zabludoff (2004) 
and Lin & Mohr (2004). 

On the theoretical front Napolitano et al. (2003) have 
used an N-body dark matter only fully cosmological sim- 
ulation of the formation of a Virgo-like cluster to make 
predictions about intracluster stars. They find that unre- 
laxed velocity distributions and (bulk) streaming motions 
of the IC stars should be common due to the large dy- 
namical timescales in clusters. Hence the evolution towards 
a relaxed coarse-grained phase-space distribution function 
through phase-mixing has proceeded much less far^ than, 
e.g., in the Galactic halo, where dynamical timescales are 
much shorter. 

The dark matter only simulations are complemented 
by various N-body simulations which invoke a more realistic 
modelling of the stellar properties of galaxies in order to 
study the fate of galaxies in a cluster environment (e.g., 
"galaxy harassment", Moore et al. 1996, "tidal stirring", 
Mayer et al. 2001 and the formation of the central cluster 
galaxy, Dubinski 1998). The advantage of this approach is 
that a very high resolution in the stellar component can 
be obtained. In relation to the properties of the IC stars, 
recent progress has been made by Dubinski et al. (2003), 
and Feldmeier et al. (2004a). 

In general the properties of the system of IC stars are 
set by two main effects: a) the cool-out of gas and subse- 
quent conversion of cold, high-density gas to stars in individ- 
ual galaxies and b) the stripping/disruption of the galaxies 
through interactions with other galaxies and the main clus- 
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ter potential. Since such interactions will generally affect the 
star-formation rate (as long as a reservoir of gas is available) 
the former process is intimately coupled to the latter. 

Only fairly recently has it been possible to carry out 
fully cosmological gas-dynamical/N-body simulations of the 
formation and evolution of galgixy clusters at a sufficient 
level of numerical resolution and physical sophistication that 
the cool-out of gas, star-formation, chemical evolution and 
gas inflows and outflows related to individual cluster galax- 
ies can be modelled to, at least some, degree of realism (e.g., 
Valdarnini 2003; Tornatore et al. 2004). Recently Murante et 
al. (2004) analyzed 117 clusters formed in a large cosmolog- 
ical TreeSPH simulation (Borgani ct al. 2004) to study the 
properties of the IC stars. They determined average density 
profiles of the IC stars as well as stars in galaxies for their 
large sample of clusters. One of the interesting conclusions 
reached is that only at fairly large cluster-centric distances 
(r > 0.4-0.5rvir) does the (spherically averaged) density of 
stars in galaxies (excluding the cD) exceed that of IC stars. 
Complementary work has recently been performed by Will- 
man et al. (2004), who simulated a Coma- like cluster at very 
high numerical resolution. They find an IC star fraction of 
10-20%, in good agreement with observations, and a ten- 
dency for the fraction to increase with time. In relation to 
on-going observational efforts, they moreover find that line- 
of-sight velocities of intra-cluster planetary nebulae can be 
used for fairly accurate cluster mass estimation, even for dy- 
namically quite unrelaxed clusters, if more than 5 fields at 
a range of projected radii are covered. 

We have recently undertaken similar simulations of 
galaxy groups and clusters. We are building on our TreeSPH 
code used for simulating galaxy formation (e.g., Sommer- 
Larsen, Gotz & Portinari 2003), improved to include mod- 
elling of non-instantaneous chemical evolution (Lia, Porti- 
nari & Carraro 2002a), metallicity-dependent, atomic radia- 
tive cooling, strong supernova, and (optionally) AGN, driven 
galactic winds and thermal conduction. The main results of 
the simulations concerning the properties of the ICM and 
cluster scaling relations are described in Paper I (Romeo et 
al. 2004) , while the characteristics of the population of the 
simulated cluster galaxies are presented in Paper II (Romeo, 
Portinari & Sommer-Larsen 2004) — in this Paper III of the 
series, we discuss our first results on the properties of the 
IC stars in two simulated clusters, with masses and temper- 
atures similar to the Virgo and Coma clusters, respectively. 
Our results complement those of Murante et al. (2004) and 
Willman et al. (2004), since we present for the IC stars 
also results on self-consistently calculated multiband sur- 
face brightness profiles, colours, abundances and detailed 
kinematics (Willman et al. 2004 also present a kinematics 
study, which we compare to in this paper). 

In section 2 we briefly describe the code and the nu- 
merical simulations, in section 3 we present and discuss the 
results obtained, and finally section 4 constitutes our con- 
clusions. 



2 THE CODE AND SIMULATIONS 

The code used for the simulations is a significantly im- 
proved version of the TreeSPH code we have used previously 
for galaxy formation simulations (Sommer-Larsen, Gotz & 
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Portinari 2003): Full details on the code are given in Pa- 
per I, here we recall the main improvements over the pre- 
vious version. (1) In lower resolution regions (which will al- 
ways be present in cosmological CDM simulations) an im- 
provement in the rmmcrical accuracy of the integration of 
the basic equations is obtained by solving the entropy equa- 
tion rather than the thermal energy equation — we have 
adopted the "conservative" entropy equation solving scheme 
suggested by Springcl & Hernquist (2002). (2) Cold high- 
density gas is turned into stars in a probabilistic way as 
described in Sommer-Larsen, Gotz & Portinari (2003). In a 
star-formation event an SPH particle is converted fully into 
a star particle. Non- instantaneous recycling of gas and heavy 
elements is described through probabilistic "decay" of star 
particles back to SPH particles as discussed by Lia, Porti- 
nari & Carraro (2002a). In a decay event a star paxticle is 
converted fully into an SPH particle. (3) Non-instantaneous 
chemical evolution tracing 10 elements (H, He, C, N, O, 
Mg, Si, S, Ca and Fe) has been incorporated in the code 
following Lia et al. (2002a,b); the algorithm includes super- 
novae of type II and type la, and mass loss from stars of 
all masses. Metal diffusion in Lia et al. was included with a 
diffusion coefficient k derived from models of the expansion 
of individual supernova remnants. A much more important 
effect in the present simulations, however, is the redistri- 
bution of metals (and gas) by means of star-burst driven 
"galactic super-winds" (see point 5). This is handled self- 
consistently by the code, so we sot k—Q in the present simu- 
lations. (4) Atomic radiative cooling depending both on the 
metal abundance of the gas and on the meta-galactic UV 
field, modelled after Haardt & Madau (1996) is invoked. We 
also include a simplified treatment of radiative transfer, by 
switching off the UV field where the gas becomes optically 
think to Lyman limit photons on scales of ^ 1 kpc. (5) 
Star-burst driven, galactic super-winds are incorporated in 
the sirrmlations. This is required to expel metals from the 
galaxies and get the abundance of the ICM to the observed 
level of about 1/3 solar in iron. A burst of star formation is 
modelled in the same way as the "early bursts" of Sommer- 
Larsen, Gotz & Portinari (2003), i.e. by halting cooling in 
the surrounding gas particles, to mimic the initial heating 
and subsequent adiabatic expansion phase of the supershell 
powered by the star-burst; this scheme ensures effective en- 
ergy coupling and feedback between the bursting star parti- 
cle and the surrounding gas. The strength of the super-winds 
is modelled through a free parameter /wind which determines 
how large a fraction of the new-born stars partake in such 
bursting, super-wind driving star formation. We find that 
in order to get an iron abundance in the ICM comparable 
to observations, /wind ^ 0.5 and at the same time a fairly 
top-heavy Initial Mass Function (IMF) has to be used. (6) 
Thermal conduction was implemented in the code following 
Cleary & Monaghan (1999). 

We ran simulations of two small and two larger galaxy 
groups as well as a Virgo-like and a (mini) Coma-like cluster 
— we shall denote the latter two clusters CI and C2 in the 
following. All 6 systems were chosen to be fairly relaoced (no 
> 1:2 merging at 1). Virial masses at z=0 were approx- 
imately 3.4x10", 7.9x10", 2.0x10" and 8.4xl0" ft^^Mo 
and X-ray emission weighted temperatures 1.1, 1.6, 3.0 and 
6.0 keV, respectively. The groups and clusters were selected 
from a cosmological, DM-only simulation of a flat ACDM 



model, with Qm~0.3, ni,=0.036, h=0.7 and (78=0.9 and a 
box-length of 150 h~^Mpc. Mass and force resolution was 
increased in Lagrangian regions enclosing the groups and 
clusters, such that mga.s=mt= 2.5x10* and mDM=1.8xl0® 
h~^MQ for the high resolution gas, star and dark mat- 
ter particles. Particle numbers are in the range 100000- 
500000 baryonic-|-DM particles. Gravitational (spline) soft- 
ening lengths of egas=e*=2.8 and eDM=5.4 /i~^kpc, respec- 
tively, were adopted. The gravity softening lengths were 
fixed ill physical coordinates from z=6 to z=0 and in co- 
moving coordinates at earlier times. To test for numerical 
resolution effects one simulation was run with eight times 
higher mass resolution and two times higher force resolution, 
yielding mgas=m,= 3.1x10^ and mDM=2.3xlO* ft^^M© and 
gravity softening lengths of 1.4 and 2.7 /i~^kpc, respectively. 

For the simulation presented in this paper /wind=0.8, 
and an Arimoto-Yoshii IMF (of slope x=-l, shallower than 
the Salpotcr slope a;=-1.35) with mass limits [0.1-100] Mq 
was adopted. AGN driven winds were not invoked. More- 
over, the thermal conductivity was set to zero assuming 
that thermal conduction in the ICM is highly suppressed 
by magnetic fields (e.g., Ettori & Fabian 2000); notice how- 
ever that, in relation to the cluster galaxy populations and 
hence also IC stars, no significant difference is found between 
simulations with zero thermal conductivity and simulations 
invoking thermal conduction at 1 /3 of the Spitzer level (Pam- 
per II). 



3 RESULTS AND DISCUSSION 

In this section we present and discuss results for clusters CI 
and C2 run at our "standard" resolution. The resolution test 
discussed in Section 3.7 indicates that none of the following 
results concerning mass and luminosity fractions, mean ages, 
surface brightness profiles, colours, abundances, kinematics 
and density profiles of the IC stars change in any significant 
way with increased numerical resolution. 

All results presented in the following refer to z=0, un- 
less another epoch is explicitly specified. At z=0 1.2x10* 
and 2.3x10* star particles are located inside of rvir=1.3 and 
2.0 h^^Mpc of the clusters CI and C2^. The corresponding 
mass in stars is 2.9x10'^'^ and 5.6x10^^ h~^MQ or a fraction 
of about 1% of the total virial mass (the stellar mass frac- 
tion is smaller for C2 than for CI - we shall return to this 
below). In both clusters the stellar mass constitutes 14% of 
the total baryonic mass inside of the virial radius. In the 
following we shall refer to star particles simply as "stars". 

When preparing the initial conditions for C2, due to 
computational limitations we sampled with split DM-I-SPH 
particles a region somewhat smaller than the "virial volume" . 
The volume outside this region, but still inside of the virial 
volume, is represented by high resolution, DM-only particles 
(with a mass scaled by 1/(1 — /6) with respect to the DM 
particles in the inner DM-I-SPH region, since they are not 

^ In this paper the virial radius rvir is the radius corresponding 
to, at z=0, an over-density 337 times the mean density of the 
Universe, as appropriate for a top-hat collapse in the adopted 
ACDM cosmology (e.g., Bryan & Norman 1998). The virial mass 
Mvir of the cluster is the mass enclosed within the region of radius 

^vir- 
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split into baryonic and DM component; is the universal 
baryon fraction, assumed to be 0.12 in this work); finally, 
the region outside the virial volume is, as usual, sampled by 
lower resolution dark matter particles with masses increas- 
ing with distance from the virial volume. As a consequence 
of this baryonic under-sampling, about half of the mass in- 
side of the virial radius of the final cluster is in the form 
of (high-resolution) DM-only particles, i.e. not split into 
SPH-I-DM. The relative percentage of these DM-only parti- 
cles increases with radius, hence the amount of late infall of 
gas and galaxies is reduced compared to the fully resolved 
case. This explains why the stellar mass fraction in C2 is 
smaller than that in CI as mentioned above (the C2 star 
mass is only about twice the CI star mass, while the respec- 
tive virial masses differ by a factor of about four). The ratio 
of the respective number of galaxies A'gai(C2)/A'gai(Cl) is 
also smaller than the ratio of the corresponding virial masses 
(again, a factor of ~2 vs. a factor of ~4, see below) and the 
average ages of the IC stars and galaxies in C2 are slightly 
larger than that in CI, due to the under-sampling of late in- 
fall of galaxies at the outskirts (Section 3.3). Since the results 
we discuss refer to the inner 1-2 Mpc, the baryonic under- 
sampling of the outskirts is not a crucial problem. Moreover, 
we have initiated a simulation of C2 with full baryonic sam- 
pling of the virial volume; at a redshift of z = 0.8 this shows 
no significant diflference in relation to the results on IC stars 
presented here. 

The stellar contents of both clusters are characterized 
by a massive, central dominant (cD) galaxy surrounded by 
galaxies and intracluster stars. The effective radii, calcu- 
lated from the B-band surface brightness profiles of the two 
cD galaxies, are Reff —5 kpc, taking the inner 80 kpc of 
the cluster as the extent of the cDs — see below. After cor- 
rection for an excess of central, young stars (section 3.2) 
the absolute magnitudes of the cDs are Mb ~ -23. Extrap- 
olating the observational "Kormendy relation" (Kormcndy 
1977) such bright cDs have Reff ~ 20h^^ kpc. Part of this 
discrepancy is due to an excess of central, young stars, as we 
shall discuss in section 3.2. We also note that in a study of 
the formation of the brightest cluster galaxy (BCG) based 
on a pure N-body simulation (but of higher resolution than 
ours) Dubinski (1998) found an effective radius of -^20 kpc 
for the central BCG. 

3.1 Identifying Cluster GalELxies and Intracluster 
Stars 

The first step in identifying the IC stars is to find all galax- 
ies in the clusters. To this end we proceed in the following 
simple way: by visual inspection of the z—0 frames the stars 
in all galaxies except the cD are located within 10-15 kpc 
from the centers of the galaxies. Guided by this typical size 
of the individual cluster galaxies, a cubic grid of cube-length 
Al=10 kpc is overlaid the cluster, and all cubes containing at 
least Nih=2 stars are identified. Subsequently, each selected 
cube is embedded in a larger cube of cube-length 3AZ. If this 
larger cube contains at least A'inin=7 star particles, which 
are gravitationally bound by its content of gas, stars and 
dark matter the system is identified as a potential galaxy. 
Since the method can return several, almost identical ver- 
sions of the same galaxy only the one containing the largest 
number of star particles is kept and classified as a galaxy. 



We tested the galaxy identification algorithm by varying Al, 
Nth, Nmin, and also the numerical resolution of the simula- 
tions (section 3.7) and found it to be adequately robust for 
the purposes of this paper. 

The algorithm identifies A'gai=42 and 94 galaxies inside 
of the virial radii of CI and C2, respectively, plus a cD at 
the center of each cluster. The ratio iVgai(C2)/iVgai(Cl) is 
smaller than the ratio of the corresponding virial masses due 
to the baryonic under-sampling in cluster C2, as discussed 
above. Apart from the cD there are no galaxies inside of 
rcD=80 kpc in either of the z=0 frames (this is not true 
in all z ~0 frames for clusters CI and C2, though, but we 
tested that the resxilts we present in the following did not 
depend on which z ~0 frame was chosen to represent the 
present epoch — see also below). 

We define the system of IC stars as the stars located at 
cluster-centric distances Tcd < r < Tvir and not inside of the 
tidal radius of any galaxy in the cluster. The tidal radius for 
each galaxy is taken to be the Jacobi limit 

where m is the mass of stars, gas and dark matter in the 
galaxy (inside of rj), D is the distance from the cluster 
center to the galaxy and M is the total mass of the cluster 
inside of D (e.g., Binney & Tremaine 1987). The cD itself is 
effectively just the inward continuation of the system of IC 
stars and the division between IC stars and cD is somewhat 
arbitrary (we hence below quote intracluster star fractions 
for rcD=80 as well as 40 kpc). The above definition of IC 
stars is conservative, since we calculate the tidal radii on 
the basis of the z=0 frame cluster-centric distances of the 
galaxies. For any galaxy which has been through at least 
one peri-center passage, this tidal radius should be taken as 
a firm upper limit. Moreover, some IC stars will be inside of 
the tidal radius of one of the cluster galaxies just by chance. 
As the total "tidal volume" of all the cluster galaxies is found 
to be a few tenth of a percent of the virial volume of the 
cluster this effect can be neglected, however. 

3.2 The Intracluster Star Fraction 

Using the above definition we find that the intracluster stars 
constitute 21 and 28% of the stellar mass inside of the virial 
radius for clusters CI and C2, respectively (adopting rcD=40 
kpc these numbers increase slightly to 28 and 34%). For 
comparison with observations it is more relevant to look at 
stellar luminosities. Since ages and metallicities are available 
for all star particles, the photometric properties are straight- 
forward to calculate treating each star particle as a Single 
Stellar Population (SSP; see Paper II for details). SSP lumi- 
nosities are computed by mass-weighted integration of the 
Padova isochrones (Girardi et al. 2002), according to the 
Arimoto-Yoshii IMF. We find that the IC stars contribute 
only 9 and 11% of the cluster B-band luminosities for CI and 
C2 respectively. This is significantly less than the observa- 
tional estimate of I/B,ic/iB,tot Jj]0.2 (e.g., Arnaboldi 2004). 
The likely reason for this discrepancy is an excess of stars 
formed fairly recently at the center of the cD and dominat- 
ing the B-band luminosity: After a period of major merging 
at z ~l-2 strong, quasi-stationaxy cooling flows develop at 
the centers of the clusters despite the strong, super-nova 
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Figure 1. Mean redshift of formation of the cD + IC stars (solid 
lines) and stars in galaxies (dashed lines) at ^=0. Results for clus- 
ter CI ("Virgo-like") are shown by thin lines and for C2 ("mini- 
Coma") by heavy lines. Statistical uncertainties are Azj <0.05. 

driven energy feedback to the IGM/ICM through galactic 
super- winds and the use of a fairly top-heavy IMF. As a 
consequence, stars continue to form steadily at the centers 
of the cDs (r <10 kpc) at rates of -350 and ~600 Mg/yr 
at z—Q for CI and C2, respectively. At z=0 the cDs are too 
blue compared to observed cDs with central B-R colours 
of ~1.0 (see Fig. 4) rather than 1.4-1.5 (note though that 
some star-formation is observed in many cDs at the base 
of the cooling flow — e.g., McNamara 2004). To make a 
(albeit crude) correction for this excess of young stars at 
the centers of the simulated cDs we temporarily discard the 
luminosity contribution of all stars in the inner rc=10 kpc 
with formation redshifts zf < Zc, where Zc is determined by 
requiring that the mean formation redshift of the remaining 
stars inside of rc should be similar to the mean formation 
redshift of stars at r>rc. We find that this is the case for 
2c ~1 and shall in the following in various cases present re- 
sults with and without this crude correction, adopting Zc=l 
for the former. With this correction I/B,ic/i'B,tot=0.29 and 
0.42 for the two clusters in better agreement with obser- 
vations (for rcD=40 kpc we obtain LB,ic/i'B,tot=0.42 and 
0.50, respectively), though somewhat high compared to the 
values found by Feldmeier et al. (2004a) for non-cD clus- 
ters. Moreover, the effective radii of the cDs increase to 9 
and 11 kpc, respectively. These values are in better agree- 
ment with observations (Kormendy 1977) and the finding 
of Dubinski (1998), but still somewhat smaller. Part of this 
remaining discrepancy may be due to the effect of gravita- 
tional "pinching" of our cDs by the excess, young central 
stars. 

Regarding the intracluster star fraction we show in Pa- 
per II that clusters CI and C2 are deficient in bright galaxies 
apart from the cDs. As discussed in Paper II, this is probably 
related to the fact that CI and C2 are fairly relaxed systems. 



for which the last major merger took place at z > 1. Hence it 
is possible that the galaxies in these clusters suffered more 
tidal stripping and merging than in an average cluster. To 
be conservative the above IC star fractions should then be 
considered upper limits to what would be found for a sta- 
tistically representative sample of cluster merging histories. 
Finally, we note that the observed IC star fraction of —20% 
has also been reproduced in the simulations of Murante et 
al. (2004) and Willman et al. (2004). 

3.3 Mean Formation Redshifts, Surface 

Brightness Profiles and Colours of the 
Intracluster Stars 

In Figure 1 we show the mean (spherically averaged) red- 
shift of formation, Zf, of the cD -I- IC stars (solid lines) and 
stars in galaxies except the cD (dashed lines) as a function 
of radial distance from the center of the cD for clusters CI 
and C2, respectively. For both clusters the average forma- 
tion redshift of the IC stars is ~3 at r > 100-200 kpc. 
The stars in galaxies (except the cD) are on average some- 
what younger with Zf^gai —2.5. This is to be expected, since 
the bulk of the IC stars originate in (proto) galaxies, which 
have been partly or fully disrupted through tidal stripping in 
the main cluster potential or by galaxy-galaxy interactions. 
In contrast, the galaxies still remaining at z=0 have poten- 
tially been able to continue forming stars out of remaining 
cold gas or gas recycled by evolved stars and subsequently 
cooled to star-forming temperatures and densities. Still, due 
to ram-pressure stripping of the hot and dilute gas reservoir 
in galactic haloes and other effects, the star-formation rate 
in the galaxies decreases significantly from z=2 to 0, consid- 
erably more so than in field galaxies cf. Paper II. Murante 
et al. (2004) find a similar trend that the IC stars are on av- 
erage older than the stars still in cluster galaxies, but they 
find somewhat lower mean redshifts of formation, Zf=1.9 
and 1.7 respectively. The reason for this discrepancy is not 
clear: it may be related to different galactic super-wind pre- 
scriptions or that we incorporate metal-dependent radiative 
cooling in the simulations, whereas Borgani et al. (2004) 
use a primordial cooling function, but is probably not due 
to numerical resolution differences (Murante 2004, private 
communication) . 

From Figure 1 it can be seen that both IC stars and 
galactic stars in cluster C2 are on average somewhat older 
than those in CI. The reason for this is numerical, rather 
than physical, related to the under-sampling of baryons in 
the outskirts of the virial volume of the cluster discussed at 
the beginning of this section: the late infall of gas and galax- 
ies onto the cluster is reduced with respect to a simulation 
with full baryonic sampling of also the outer parts of the 
virial volume. 

In Figure 2 we show for cluster C2 the azimuthally av- 
eraged UBVRIJK surface brightness profiles of the cD-|-IC 
stars, both with and without the correction for central, 
young cD stars discussed above (the H-band results are 
omitted for clarity). The orientation is such that the cD 
is viewed along the minor axis, but other orientations give 
very similar results. The light profiles are approximately de- 
scribed by r^''* laws. In reality the slope flattens with in- 
creasing r, so the system can be seen as a central dominant 
elliptical galaxy surrounded by an extended envelope, with 
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Figure 2. Multiband (UBVRIJK going bottom up) surface 
brightness profiles of cD + IC stars for cluster C2, shown with 
(thin lines) and without (thicker lines) the correction for the ex- 
cess young, central stars, discussed in the text. Also shown (thick 
dotted line) is V-band surface photometry for the cD in the rich 
cluster A1413 (T ~8 keV) obtained by Feldmeier et al. 2002 



the slope of the surface brightness profile becoming close to 
constant beyond r ~40 kpc. These are characteristics of ob- 
served cD galaxies (e.g., Feldmeier et al. 2002). Observation- 
ally, cD-|-IC stars can be traced by surface photometry to 
Yiim —28.3 mag/arcsec'^ (e.g., Feldmeier et al. 2002), which 
for cluster 02 corresponds to r ~350 kpc. In Figure 2 we 
also show V band surface photometry of the cD -|- envelope 
in the rich cluster A1413 (T ~8 keV) obtained by Feldmeier 
et al. (2002). Though the observational V band profile is 
slightly fiatter than what we predict for the V band, overall 
agreement is quite good, especially when it is taken into ac- 
count that the sampling of the baryons becomes increasingly 
incomplete with radius in cluster 02 (see above). In Figure 
3 we show the surface brightness profiles of 01 ("Virgo"). 
Also shown is two B-band data points derived for Virgo 10 
stars (indirectly) using planetary nebulae by Arnaboldi et 
al. (2002). Our predicted B-band profile fits the inner point 
well, but not the outer one, which is too high (in fact slightly 
higher than the inner one). The Virgo cluster is known, how- 
ever, to have an unusually fiat galaxy surface brightness pro- 
file and possibly be dynamically unrelaxed (Binggeli et al. 
1987). 

Shown in Figure 4 are the azimuthally averaged B-V, 
B-R and I-K colours of the cD-|-IO stars, again with and 
without the correction discussed above. At about 10 kpc 
the colours (e.g., B — V ~0. 90-0. 95) are typical of the stellar 
populations in the inner parts of cDs (e.g., Mackie 1992). 
Within 10 kpc this is also the case when correcting for the 
excess of central, young stars — if not, the core of the cD 
gets too blue with central colours approaching those of spi- 
ral galaxies. Outside of 10 kpc the colour of the cD-|-IO stars 
gets bluer with increasing R, but the gradient is very shal- 



Figure 3. Same as in Figure 2, but now for cluster CI. Data 
points are B-band IC star estimates based on planetary nebulae 
observations in the Virgo cluster by Arnaboldi et al. 2002 
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Figure 4. Azimuthally averaged B-V (solid lines), B-R (dashed 
lines) and I-K (dot-dashed lines) colours of cD -|- IC stars for 
clusters C2 (thick lines) and CI (thinner lines), respectively. The 
results of applying the correction for the excess, central young 
stars is shown with thin dotted lines. 

low: A(B-V) per kpc of ~ -0.0003 As the mean age of the 
stars is approximately constant with R (see Fig. 1), this is 
mainly due to a decrease in metal abundance with R (see 
below). At 100 kpc B-R ~1.45, typical of sub-L* E and SO 
galaxies (e.g., Gladders et al. 1998). We shall discuss this in 
a forthcoming paper, but note that cDs are found in general 
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Figure 5. Spherically averaged iron abundance of cD+IC stars 
(solid lines) and stars in cluster galaxies (dashed lines) for clusters 
CI (thin lines) and C2 (thick lines), respectively (results for r <20 
kpc are not shown for clarity). 



to have quite flat colour profiles, in some cases getting red- 
der with R (e.g., Mackie 1992; GariUi et al. 1997; Gonzalez 
et al. 2000). 



3.4 Abundance Properties of the Intracluster 
Stars and Cluster Galaxies 

In Figure 5 we show the spherically averaged iron abundance 
of the cD+IC stars as well as of the stars in cluster galax- 
ies as a function of cluster-centric distance. Iron is super- 
solar in the cD-|-IC stars at r < 100 kpc and decreases to 
about half solar at large r; the stars in cluster galaxies fol- 
low a similar trend. The fairly large overall iron abundance 
of the IC stars, as compared to, e.g., stars in the halo of the 
Milky Way, reflects that it is the galaxies (past and present) 
which have to enrich the 5-10 times more massive hot in- 
tracluster medium (ICM) to an iron abundance of about 
1/3 solar (Paper I). DurreU et al. (2002) carried out HST 
observations of an IC field in the Virgo cluster at an aver- 
age projected distance of 150 kpc from M87 (which for the 
purposes here can be assumed coincident with cluster cen- 
ter). They confirm an excess of red number counts, which 
they interpret as IC RGB stars. By comparison with obser- 
vations of a dwarf irregular, they conclude that these stars 
have -0.8<[Fe/H]<-0.2. Though this is somewhat less than 
we predict at r ~150 kpc, the discrepancy is only a factor 
of about two, and the observational abundances are clearly 
significantly larger than those of the halo stars of the Milky 
Way. 

Figure 6 shows the corresponding oxygen-to-iron ratios 
as a function of r. [O/Fe] is super-solar everywhere. This is 
in agreement with present estimates for the luminous ellip- 
tical galaxies that contain most of the stellar mass in clus- 
ter galaxies. For the IC stars no observational information 



Figure 6. Spherically averaged oxygen-to-iron abundance ratio 
of cD+IC stars (solid lines) and stars in cluster galaxies (dashed 
lines) for clusters CI (thin lines) and C2 (thick lines), respectively. 
The limiting case for pure SNII enrichment [O/Fe]=0.55 is shown 
by the thin, dotted line (results for r <20 kpc are not shown for 
clarity) . 



is currently available (Arnaboldi 2004, private communica- 
tion). 

For pure type II supernovae enrichment and with the 
Arimoto-Yoshii IMF, one expects [O/Fe]sNii=0.55 (e.g., 
Lia, Portinari & Carraro 2002a), so it follows from Fig. 6 
that SNe la do contribute somewhat to the enrichment of 
the cD+IC stellar populations, and (not surprisingly) even 
more so for the stars still in galaxies at z=0 (in fact for 
an Arimoto-Yoshii IMF the global {t oo) value of the 
SNII-l-SNIa enrichment is [O/Fe]=0.18). 



3.5 Kinematics of the Intracluster Stars and 
Cluster Galaxies 

Using observed velocities of planetary nebulae it will ulti- 
mately be possible to kinematically "dissect" the systems of 
cD-|-IC stars in nearby galaxy clusters, such as Virgo. It is 
hence of considerable interest to determine for our simula- 
tions the velocity distribution of the cD-|-IC stars and com- 
pare it to that of cluster galax;ies. To this end we proceed 
as follows: At z—0 cluster CI and the cD at the center are 
somewhat flattened (ellipticity e = 1 — b/a ~0.4) with sim- 
ilar minor axis orientations (we shall denote the minor axis 
of the cD the "z-axis" in the following). Cluster C2 is only 
slightly flattened (e —0.2), and so is the cD at the center, 
with approximately perpendicular minor axis orientations. 
For each cD-|-IC star and each cluster galaxy we determine 
three perpendicular velocity components: The radial compo- 
nent Vr=v ■ Cr, where eV is the unit vector pointing radially 
away from the center of the cluster, the perpendicular (tan- 
gential) component v^=v ■ e^, where is the unit vector 
perpendicular to eV and aligned with the x-y plane and the 
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Figure 7. For cD+IC stars in cluster CI we show tlie velocity 
dispersions ar (solid lines), cr^ (dotted line) and erg (dashed line). 
The statistical uncertainties are marked with small errorbars. The 
similar quantities are shown for the cluster galaxies by filled cir- 
cles, crosses and filled triangles, with larger errorbars marking the 
statistical uncertainties. 
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Figure 8. Same as Figure 7, but for cluster C2. 

third (tangential) component vg—v- eg, where eg is the unit 
vector eg = eV X e^t. 

We calculate the mean rotation and velocity disper- 
sions (Tr, and ae of cD-|-IC stars and galaxies in spherical 
shells. Outside r ~10 kpc rotation is found to be dynam- 
ically insignificant with i)0< 20-40 km/s in both clusters. 
When correcting for the excess young, central stars rota- 



Figure 9. Velocity anisotropy parameter 7 for cD+IC stars (solid 
line) and galaxies (solid circles) for cluster CI. The same is shown 
for cluster C2 by dashed line and solid triangles. An isotropic 
velocity distribution has 7=1 indicated by a thin, dotted line. 

tion is found to be also dynamically unimportant at r <10 
kpc. In Figures 7 and 8 we show the velocity dispersions 
of the cD-l-IC stars and cluster galaxies in clusters CI and 
C2, respectively, versus cluster-centric distance. The cD-|-IC 
stars are kinematically significantly colder than the cluster 
galaxies: The velocity dispersions of the cD-f IC stars are in 
the inner parts of the clusters (r ~100-500 kpc) only about 
half of those of the cluster galaxies, increasing slightly to 
about 70% at r ~l-2 Mpc. As stars and galaxies are moving 
in the same gravitational potential this implies by Jean's 
theorem (e.g., Binney & Tremaine 1987) that the (number) 
density distribution of cluster galaxies is significantly fiatter 
than that of cD-fIC stars, in particular in the inner parts 
of the clusters — see below. Willman et al. (2004) also find 
that the IC stars are kinematically colder than the cluster 
galaxies, except for the innermost part of the cluster. The 
difference here is probably due to the different definitions of 
IC stars used: Willman et al. (2004) require the IC stars to 
be unbound from the "cD" also — this results in large IC 
star velocity dispersions in the inner parts of the cluster. 

The velocity distribution of the cD-|-IC stars in cluster 
CI is highly radially anisotropic at r > 100 kpc, and, within 
the statistical uncertainties, this is also the case for the clus- 
ter galaxies. We quantify this in Figure 9, where we show 
7 = •y/(2(T^/(CT^ + ag)), where 7^ is the ratio of the kinetic 
energy in radial and mean tangential (ID) motions, respec- 
tively. For an isotropic velocity distribution 7=1. In cluster 
CI there is almost twice as much kinetic energy in the radial 
direction than in the sum of the two tangential directions at 
r > 100 kpc. Interestingly, within the statistical uncertain- 
ties the same appears to hold for the galaxies, which at first 
seems surprising, since one would expect that it is predom- 
inantly galaxies on radial orbits which are tidally disrupted 
and transformed into cD-f IC stars. The explanation for this 
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Figure 10. Projected velocity dispersions for cluster CI of 
cD+IC stars (thick lines) and dark matter (thinner lines) along 
the minor axis (z-axis; dashed lines) and averaged along the x 
and y-axis (solid lines) of the cD/cluster. Thin dotted lines mark 
the effect of making the correction for the excess young central 
stars. Statistical uncertainties are A(t<20 km/s. 



lack of difference in the orbital characteristics of the two 
populations is likely the combined effect of radial infall of 
and dynamical friction on cluster galaxies, but we defer a 
detailed discussion of this to a forthcoming paper. We note 
that Willman et al. (2004) find the velocity distribution of 
the IC stars to be more radially anisotropic than that of the 
galaxies in the outer parts of their simulated cluster. In clus- 
ter C2 the velocity distribution of the cD+IC stars is only 
mildly radially anisotropic and the velocity distribution of 
the cluster galaxies is, within the statistical uncertainties, 
isotropic. We also note that Sommer-Larsen et al. (1997) 
found that the velocity distribution of Galactic halo stars 
becomes tangentially anisotropic in the outer halo — the 
differences probably relate to different accretion kinematics 
during the build-up of the various haloes — this will also be 
discussed in a forthcoming paper. 

Finally, concerning the flattening of the systems for 
cluster CI is about 50% larger than (jg for cD+IC stars 
as well as galaxies, indicating that these systems are flat- 
tened by anisotropic velocity distributions — see below. For 
cluster C2 ~ ere, consistent with this cluster being only 
mildly flattened. 

Observationally, for the cD-|-IC stars one will only be 
able to determine line-of-sight velocities using planetary 
nebulae, not full 3D velocities. For direct comparison with 
observations we show in Figure 10 the projected velocity 
dispersions of the cD-|-IC stars and (for comparison) of the 
dark matter in cluster CI versus projected distance from the 
center of the cD. By dashed lines we show the azimuthally 
averaged velocity dispersions along the minor axis of the cD. 
The projected stellar velocity dispersion is 350-400 km/s at 
the center of the cD, increases to ~500 km/s at K ~10- 



30 kpc and then decreases gradually with increasing R to 
about 300 km/s at -Rvir- An increase of the stellar velocity 
dispersion with R has been observed in some cDs, such as 
A2029 (Dressier 1979), and has been interpreted as marking 
the transition from galactic to intracluster stars. The effect 
is less pronounced after correction for the excess young cen- 
tral stars, but we note that these stars still contribute to the 
gravitational potential and therefore boost the central veloc- 
ity dispersions. The projected velocity dispersion of the dark 
matter follows a similar trend with R as that of the stars, 
but is significantly larger. As the stars and dark matter are 
moving in the same gravitational potential this implies that 
the density distribution of dark matter is significantly flatter 
than that of the cD-f IC stars — see below. 

The projected velocity dispersions along the major axis 
of the cD are shown by solid lines. They have been calculated 
as the average of the projected velocity dispersions along the 
X-axis and y-axis. Along the x-axis the velocity dispersions 
are calculated for stars and dark matter particles projecting 
onto ±30 deg. wedges along the positive and negative y di- 
rections, and vice versa for dispersions along the y-axis. The 
major axis velocity dispersions follow similar trends as the 
minor axis ones. However, for the cD-|-IC stars it is signifi- 
cantly larger and for the dark matter halo somewhat larger 
than along the minor axis. This shows that both systems are 
fiattened by anisotropic velocity dispersions and, as stated 
above, oriented in similar ways. For cluster C2 the findings 
are similar (not shown) though the cD is not aligned with 
the cluster in this case and both systems are only mildly 
fiattened. 

Summarizing our findings for the Virgo-like cluster, the 
typical projected velocity dispersion for the cD-|-IC stars 
at ii>50 kpc is 300-600 km/s depending on orientation 
and projected distance from the cluster center. Freeman et 
al. (2000) find a (projected) velocity dispersion of 752±110 
km/s for 23 IC planetary nebulae in the Virgo-cluster at a 
projected distance of about 150 kpc from M87. This value is 
somewhat larger than our prediction, but it may be due to 
the "real" Virgo-cluster being kinematically unrelaxed (see 
also section 3.3); the observational sample is too small to 
test for velocity substructure (Freeman 2004, private com- 
munication) . 

Since clusters are dynamically fairly young it is impor- 
tant to establish that the results on kinematics etc. are ro- 
bust, i.e. they do not depend on the particular frame chosen 
for the analysis. To this end we analyzed for each cluster two 
frames 1.5 and 3 Gyr before the present time, corresponding 
to redshifts 0.12 and 0.26, respectively. We find that none of 
the results presented in this paper change in any significant 
way when going to these earlier frames. 

3.6 Density Distributions of the Intracluster 

Stars, Cluster Galaxies and the Dark Matter 

In Figure 11 we show for cluster CI the spherically averaged 
density distributions of the cD-|-IC stars (with and without 
the correction for the excess, young central stars), the stars 
in cluster galaxies, the dark matter and, with an arbitrary 
normalization, the number density of cluster galaxies. As 
suggested previously, the number density profile of cluster 
galaxies in the inner parts of the cluster as well as the dark 
matter density profile in general are significantly flatter than 
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Figure 11. Spherically averaged density profiles for cluster CI of 
cD+IC stars (solid line), dark matter (dot-dashed line), stars in 
cluster galaxies (dashed line) and, with arbitrary normalization, 
number density of cluster galaxies (dotted line). The result of 
applying the correction for the excess, young central stars is shown 
by the thin, dotted line. 

the density profile of cD+IC stars. For the cluster galaxies 
this is a manifestation of the transformation of galaxies to 
cD+IC stars by tidal disruption in the inner parts of the 
cluster. 

Moreover, it follows from Fig. 11 that the average (mass) 
density of stars in cluster galaxies only exceeds that of the 
IC stars beyond about half the virial radius - this is in good 
agreement with the findings of Murante et al. (2004) . 

3.7 A Numerical Resolution Test 

It is important to test whether the properties of the IC star 
population presented in this paper depend much on the nu- 
merical resolution of the simulations. In the present paper, 
we address resolution effects by discussing a smaller sys- 
tem, a group of virial mass 3.4 x lO^^ft~^M0 and emission- 
weighted temperature 1.1 keV. The group was run at the 
same "standard" resolution as clusters CI and C2, as well as 
with eight times higher mass and two times higher force res- 
olution. This results in particle masses of mgas=m,=3.1xl0^ 
and mDM=2.3xlO* /I'^M© and gravity softening lengths of 
1.4, 1.4 and 2.7 /i~^kpc, respectively. In order to enable an 
optimal comparison between the normal and high resolution 
runs only Fourier modes up the Nyquist wavenumber of the 
normal resolution simulation were used to prepare the ini- 
tial conditions for the high resolution run (i.e., additional 
high-wavenumber modes up to the Nyquist wavenumber of 
the high resolution simulation were not added to the Fourier 
modes) . 

The resulting group looks like a "scaled down" version 
of clusters CI and C2. In particular, it has a prominent 
central galaxy, which we denote as the "cD", like in the 



Figure 12. For a A/vir=3.4xl0^^ h'^MQ (T=l.l keV) group is 
shown at 2=0.1 the cumulative mass in cD+IC stars (solid lines) 
and in stars in cluster galaxies (dashed lines) outside of r=30 kpc 
for the normal resolution (thin lines) and high resolution (thick 
lines) simulations, respectively. 



cluster simulations — in this respect it is similar to the so- 
called "fossil" groups (e.g., Jones et al. 2003). 

In Figure 12 we show the cumulative mass of cD-f IC 
stars and stars in galaxies outside of r=30 kpc for the nor- 
mal and high resolution runs at z—Q.l (at z—O one large 
group galaxy is merging with the cD). There is reasonable 
agreement between the runs — the somewhat larger mass in 
IC stars in the high resolution simulation at r > 200-300 kpc 
is likely due to the better resolution of star-forming gas in 
the lower over-density regions which come to populate the 
outer parts of the group. 

As one goes to higher resolution one might naively ex- 
pect that the fraction of IC stars in the simulations decreases 
as the cluster galaxies are increasingly well resolved. The in- 
teresting indication from Figure 12 is, however, that if any- 
thing the fraction of IC stars slightly increases with increas- 
ing numerical resolution. 

Shown in Figure 13 at 2=0.1 are the azimuthally aver- 
aged velocity dispersions of cD-|-IC stars and dark matter 
projected along the minor axis of the cD for the normal and 
high resolution simulations, and with the correction for ex- 
cess young central stars discussed previously applied. The 
agreement is overall reasonable, especially when it is taken 
into account that there are at z=0.1 just 1235 cD-fIC stars 
inside of the virial radius in the normal resolution simulation 
of this group. This is a factor 10-20 less than the correspond- 
ing numbers for clusters CI and C2. The increase in central 
velocity dispersion (-R^ 10 kpc) for the high resolution sim- 
ulation, in particular for the stars, is likely an effect of the 
increased force resolution. Note also the qualitative similar- 
ity between the velocity dispersion profiles shown in Fig. 13, 
and the results for cluster CI (Fig. 10). 

Finally, as a test of resolution dependence of chemical 
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Figure 13. Velocity dispersions projected along the minor axis 
of the "cD" for the cD+IC stars (lower curves and symbols) and 
dark matter (upper curves and symbols) in the group. Results for 
the normal resolution simulation are connected by dashed lines 
and are shown as solid circles with errorbars (stars) and solid 
triangles with errorbars (dark matter). Results for the high reso- 
lution simulation are shown by solid lines. 

evolution and metal abundances we show in Figure 14 the 
cumulative mass of Iron in cD+IC stars and stars in galaxies 
outside of r=30 kpc for the normal and high resolution runs. 
There is reasonable agreement between the runs and we note 
that in relation to Figure 12 this is a nontrivial result, since 
a significant fraction of the iron produced by the stars ends 
up in the hot ICM (Paper I). 

In general, as mentioned at the beginning of section 
3, we carry out all the analysis of the (normal resolution) 
cluster simulations CI and C2 presented in this paper for 
the normal and high resolution group simulations also. On 
this basis we find that all results given in this paper appear 
largely robust to resolution changes. Ultimately, however, 
this can only be properly checked by running the "produc- 
tion" simulations CI and C2 at higher numerical resolution. 
To this end we have initiated a high resolution simulation of 
CI (with similar numerical characteristics as the above high- 
resolution group simulation) — at z=1.3 we find no signifi- 
cant differences between the "standard" and high-resolution 
simulations of CI in relation to the results presented in this 
paper. 



4 CONCLUSIONS 

In this paper we have discussed the properties of the in- 
tracluster (IC) stellar populations. Our results are based 
on cosmological simulations of galaxy clusters including 
self-consistently metal-dependent atomic radiative cooling, 
star-formation, supernova driven galactic super-winds, non- 
instantaneous chemical evolution and the effects of a meta- 
galactic, redshift dependent UV field. In relation to mod- 



Figure 14. Same as Figure 12, except that cumulative masses of 
Iron are shown. 



elling the properties of the IC stars this is an important 
step forward with respect to previous theoretical works on 
the subject which were based on purely N-body (collision- 
less) simulations, apart from the recent, interesting work by 
Murante et al. (2004). 

The main results from our simulations of a Virgo-like 
(CI) and a (sub) Coma-like (C2) galaxy cluster regarding 
the IC stellar populations are as follows: 

The intracluster (IC) stars are found to contribute 20- 
40% of the total cluster B-band luminosity at z=0 and to 
form at a mean redshift Zf ~3, somewhat larger than the 
mean formation redshift of the stars in cluster galaxies which 
is about 2.5. This difference corresponds to a time span of 
about 0.5 Gyr. Murante et al. (2004) find somewhat lower 
mean formation redshifts of 1.9 and 1.7, respectively — the 
reason for this discrepancy between the two works remains 
to be identified. 

We calculate UBVRIJHK surface brightness profiles of 
the IC star populations and find that the profile of the larger 
cluster matches the observed V-band profile of Abell 1413 
(r ~8 keV) quite well. For the Virgo-like cluster we fail, 
however, to match the flat profile between 200 and 400 
kpc projected cluster-centric distance observationally (but 
indirectly) inferred for the Virgo cluster by Arnaboldi et 
al. (2002) from planetary nebulae counts. We note though 
that the galaxy distribution in Virgo is also unusually flat 
(Binggeli et al. 1987). The V-band surface brightness is 
found to reach Yum —28.3 mag/arcsec^ (Feldmeier et al. 
2002) in cluster C2 at r ~350 kpc and in cluster CI at 
r ~250 kpc. The typical colour of the IC stellar population 
is B-R=1.4-1.5, comparable to the colour of sub-L* E and SO 
galaxies. The mean Iron abundance of the IC stars is about 
solar in the central part of the cluster (r ~100 kpc) gently 
decreasing to about half solar at the virial radius. The IC 
stars are a-element enhanced with, e.g, [O/Fe] increasing 
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slightly with r and characterized by a typical [O/Fe]~0.4 
dex. 

The IC stars are kinematically significantly colder than 
the cluster galaxies: The velocity dispersions of the IC stars 
are in the inner parts of the cluster (r ^^100-500 kpc) only 
about half of those of the cluster galaxies increasing slightly 
to about 70% at r=l-2 Mpc. The typical projected velocity 
dispersion in the Virgo-like cluster at -R > 50 kpc is 300- 
600 km/s depending on orientation and projected distance 
from the cluster center. Rotation is found to be dynamically 
insignificant for the IC stars. The velocity distributions of IC 
stars and clusters galaxies are in one cluster highly radially 
anisotropic, in the other close to isotropic. 

A test simulation of a T ~1.1 keV group at higher nu- 
merical resolution indicates that the results presented are 
largely robust to resolution changes. Work is in progress to 
simulate also at higher resolution the two clusters discussed 
in this paper. Moreover, we are in the process of enlarging 
significantly our sample of galaxy clusters. 
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